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1. Introduction 

Obtaining precise physical results from lattice calculations requires a well controlled contin- 
uum limit and, for many quantities, non-perturbative renormalisation. Ideally the renormalisation 
scheme should be not only non-perturbative but also mass independent and preferably gauge in- 
variant. Schrodinger functional (SF) schemes |[TJ, ||[ |5|] are known to fulfill these properties. Addi- 
tionally, to ease the burden of taking the continuum limit, 0(a) improvement is highly desirable. 
However, to eliminate the many counterterms necessary when applying the standard 0(a) improve- 
ment program with Wilson fermions, we would like to capitalize on the automatic 0(a) improve- 
ment provided by maximally twisted mass fermions (see [|]] for a review). Unfortunately, bulk 
automatic 0(a) improvement with Wilson fermions and the standard SF (sSF) boundary conditions 
(BCs) are not compatible. 0(a) improvement is only possible introducing a number of additional 
bulk improvement counter-terms to the action and operators. Since there are extensive calculations 
with maximally twisted mass fermions [^, 0] it would be clearly desirable to employ the SF scheme 
while keeping automatic 0(a)-improvement. 

A new formulation of the SF has been developed in Ref. which we will refer to as the 
chirally rotated SF (#SF), that implements a SF scheme while maintaining automatic 0(a) im- 
provement for massless Wilson fermions. The #SF is related (in the continuum) to the sSF by 
means of a non-singlet chiral transformation, i.e. they are equivalent in the continuum limit. How- 
ever, when using massless Wilson fermions as a lattice regulator, #SF BCs are invariant under a 
subgroup of the chiral symmetry transformations broken by the Wilson term (in contrast to sSF 
BCs). As a result #SF BCs are compatible with automatic 0(a) improvement. 

The three-dimensional boundaries of the SF lead to an unavoidable dimension four boundary 
operator. Additionally, regulating the #SF with Wilson fermions induces the usual bulk mass 
operator as well as a dimension three boundary operator. The dimension four boundary operator 
is irrelevant, and hence the corresponding coefficient can be safely fixed by perturbation theory in 
order to eliminate the corresponding 0(a) boundary contributions. The bulk operator is relevant 
and is handled by the standard non-perturbative tuning of the bare quark mass, equivalently k, to 
its critical value. The dimension three operator is also relevant and can spoil not only the automatic 
0(a) improvement but also the universality of the continuum limit. This requires an additional 
non-perturbative tuning of one more counterterm, Zf. However, having tuned both k and zu all 
operators are automatically 0(a) improved and no further counterterms are necessary. 

Here we present the non-perturbative tuning of K and Zf for the #SF in the quenched approx- 
imation. We demonstrate the feasibility of tuning both parameters simultaneously. In particular, 
the inclusion of the bulk dimension five operator, with corresponding counterterm c sw , as used in 
Ref. is found to be unnecessary. 

2. Boundary conditions 

The #SF is related to the sSF by a non-singlet chiral transformation, % = ex P( — ^75^ 3 /4)i//, 
where y is the fermion doublet in the Nf = 2 standard formulation, % is the corresponding doublet 
in the rotated basis and T 3 is a Pauli matrix. This field transformation maps the sSF BCs to the ^SF 
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BCs, 



e+zWk=o=o 
zWe+k=o = o 



G-zWU=r = o 

X(x)Q-U=t = 0, 



(2.1) 



where T is the Euclidean time extent and Q± are projectors given by 

Q± = \ (l±/7o75T 3 ) • 

Thus the Q± are simply the chirally rotated projectors corresponding to the sSF projectors, P± = 
1/2(1 ± 70). However, once the theory is regularised on the lattice, we must ensure that the BCs 
in ( |Q| ) are in fact recovered in the continuum limit. Using orbifolding techniques, it was shown 
that the BCs can be implemented at finite lattice spacing by a simple modification of the standard 
Wilson-Dirac operator, Dw, near the time boundaries [[H]]. The resulting action is 



5 = fl 4 XLlW(% + 'no)z(x) 

,v =0 x 



(2.2) 



and the modified Wilson-Dirac operator is given by 
where K is the time-diagonal contribution to Dw 



' -U(x,0)P-X(x + ad) + (aK + iY 5 T 3 P_)x(x) if x = 

aD w xi x ) if 0<x <r 

(aK + iY 5 z 3 P + )x(x)-U t (x-a6,0)P + x{x-a6) if x = T 



(2.3) 



3. Boundary counterterms 

To ensure the correct continuum limit, we must account for all relevant operators allowed by 
the symmetries of the action above. This means dimension four or less for the bulk action. There 
is one such operator, ~%%, and the corresponding counterterm is the term proportional to the critical 
quark mass, m cr , or equivalently K^ r . This is the standard operator that is present for all Wilson 
actions due to the breaking of chiral symmetry by the Wilson term. 

Similarly, we must include all permitted boundary operators of dimension three or less. Again, 
the one allowed operator is XX which gives rise to the following counterterm to the lattice 
action, 

8S 3 = (z f -l)a 3 £ {XX\xo=o+XX\xo=t) ■ 

X 

Such an operator would be forbidden in the continuum action, but the reduced symmetries of 
the Wilson action do not allow us to exclude this operator on the lattice. The presence of 553 
can then be understood as necessary to restore the symmetries broken by the Wilson term in the 
continuum limit. The fact that it is a relevant operator implies that we must compute the bare 
coupling dependence of Zf non-perturbatively, just as for K. 
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Furthermore, we must examine those irrelevant operators that lead to 0(a) contributions. In 
the bulk, there is the dimension five Sheikholeslami-Wohlert term, but automatic 0(a) improve- 
ment eliminates the need for this operator. Yet, there does remain an 0(a) contribution from the 
boundary due to the irrelevant dimension four operator 

8S 4 = (d s - l)a 4 £ (XYkDa\x =o + XYkD k z\ Xo =T) ■ 

X 

Such a contribution is present in all SF formulations [J3j] and is not due to the particular lattice 
action or BCs we have chosen. In fact, d s plays a role that is analogous to the c, counterterm in the 



sSF [11]. Given that 8S4 is an irrelevant operator, d s can be computed in perturbation theory. For 



the investigation presented here, we simply use the tree-level value of 1 /2. 

4. Tuning conditions 

The non-perturbative determination of K and Zf requires imposing conditions at finite lattice 
spacing that ensure the restoration of all expected symmetries in the continuum limit: parity and 
flavour symmetries in the % — basis 1 . Moreover, these conditions should be imposed at each lattice 
spacing while fixing a suitable renormalised quantity. In this work, we keep the renormalised SF 
coupling, g, fixed. This is equivalent to fixing the physical size of the box, L. All other dimensionful 
quantities must scale with L, so we choose T = L, evaluate all correlation functions at xq = T /2 
and use periodic boundary conditions with 8=0. 

Before specifying the tuning conditions, we define the following boundary to bulk correlation 
functions 

<(*o) = -(K(x)£i) gf ± (x Q ) = -(P a (x)£ b ± ) 
where the boundary operator, is defined for the xq = boundary by 

^ = a^X(y)Y5^ a Q±^)e rp( "- t] , 

y,z 

the bulk operators A a ^(x) and P a (x) are the axial current and pseudoscalar density in the ^-basis, 
and the boundary fields for xq = are defined as 

£(x) =U(x -a,x;0)x(x)\x =a =X( x )U t (x Q - a,x;0)\ Xo=a . 

To tune K to its critical value, we adopt the standard procedure of imposing a vanishing PCAC 
mass. To tune Zf, we require the y 5 Ti-odd correlation function g 1 ^ to vanish, 

_ ^d att ^A_( r / 2 ) „ _ h/^/on n ,< n 

mpcAC = - UTFT^s = <?A— = g A _ (T/ 2 ) = • ( 4 -l) 
2 §p ( T /2) 



The second condition in particular is sensitive to the symmetries broken by the lattice action (|2.2|), 
and both conditions together ensure that in the continuum limit all broken symmetries are indeed 
restored. Imposing different symmetry restoration conditions would give rise to different values of 
K and Zf that would differ amongst themselves by cutoff effects. It will be important to study the 
sensitivity of K and zf to the particular definitions used in order to better understand the intrinsic 
uncertainty in the determination of these counterterms. 

1 We recall that in the %— basis parity and flavour symmetries take a slightly different form (see ref. |^|| for a 
discussion about the dependence of the symmetries on the basis adopted). 
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L/a 


P 


4 (*sf) 


CfSF) 


Kb (sSF) 


Tuning at a hadronic scale, /I ~ 300 MeV 


8 


6.0219 


1.8090 (32) 


0.153530(24) 


0.153371(10) 


10 


6.1628 


1.7920 (30) 


0.152134(17) 


0.152012(7) 


12 


6.2885 


1.7664 (51) 


0.150815(22) 


0.150752(10) 


16 


6.4956 


1.7212(83) 


0.148945 (25) 


0.148876(13) 




Tunin 


g at an intermediate scale, /I ~ 1 GeV 


8 


7.0197 


1.5467(15) 


0.144501(13) 


0.144454(7) 


12 


7.3551 


1.5126 (23) 


0.143113(12) 


0.143113(6) 


16 


7.6101 


1.4942 (37) 


0.142112(13) 


0.142107 (6) 


Tuning at a perturbative scale, jj. ~ 30 GeV 


8 


10.3000 


1.29730(67) 


0.1354609(54) 


0.135457 (5) 


12 


10.6086 


1.2954(11) 


0.1351758(56) 


0.135160(4) 


16 


10.8910 


1.2858(15) 


0.1348440 (61) 


0.134849 (6) 



Table 1: Tuning results at a hadronic, intermediate and perturbative scale. We give the critical values, Zf 
and rc cr , calculated in this work for the ^SF. For reference, we also give JC cr for the sSF jl2| , pj[ [j4j]. 

5. Tuning results 

To check the practicality of tuning both K and Zf non-perturbatively for the #SF, we perform 
the tuning at three values of the renormalisation scale pL = 1/L, corresponding to a hadronic (g 2 
fixed with L = 1.436ro), an intermediate (g 2 = 2.4484) and a perturbative (g 2 = 0.9944) scale. The 
results at these three points are summarised in Tab. ||. We now briefly explain the procedure we 
used to perform the tuning, showing examples from our most difficult point at the hadronic scale 
and for the smallest lattice, L/a = 8. 

The values of j8 used are given in Tab. |l]and are taken from Ref. [jFJ]]. The tuning is performed 
in several steps. First, we calculate rapcAC and gA- at four values of zt, and for each value of Zf, we 
use four values of K, thus giving 16 pairs of fc and Zf. This allows us to determine gA_ as a function 
of mpcAC for each value of zt, as illustrated in Fig. [IJ. For each value of Zf, we perform a linear 
interpolation of g&- in terms of mpcAC to the point wipcac = 0. This determines the values of gA- 
at mpcAC = for each of the four values of Zf, as shown in Fig. ||. We now interpolate these values 
of gA- as a function of Zf to the point of vanishing g&_, thus giving us the critical value z,f. 

Next we determine K^ r . Using the same 16 pairs of K and Zf, we calculate wipcac as a function 
of K for each Zf. This is shown in Fig. ||. Note that rapcAC has a very mild dependence on zt, so 
the four curves at fixed Zf are nearly indistinguishable. Interpolating in fc to the point of vanishing 
PCAC mass, we obtain the critical values of k at each z,f- The resulting values of K as a function of 
Zf are shown in Fig. ||. We now interpolate these results in Zf to the previously determined value of 
Zf, thus determining the value of Kc r . 

A key observation of this work is the mild dependence of mpcAC on Zi, at least in the region 
near K^ r and Zf. You can easily see this in Fig. [| The consequence of this is clear in Fig. ^ 
the determination of K^ r also has a weak dependence on z,f and the errors of both are relatively 
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tc-2.951i.15) 



L=1.436r | , 
8 4 , B=6.0219 



Figure 1: Plot of gA- versus mpcAC- 



Figure 2: Plot of gA- versus Zf ■ 



I,. 



'"X. 



L=1.436,„ 
-0 01 - 8* 6=6.0219 



0.153 0.15325 0.1535 



"'1 - 



5^1.8090(32) 
K —0. 153530(24) 



L=1.436,„ 
8, 8=6-0219 



' 1.785 1.79 1.795 1.8 



1.805 1.81 1.815 1.82 1.825 



Figure 3: Plot of mpcAC versus K. 



Figure 4: Plot of Kc T versus zt ■ 



independent. If this behaviour persists with dynamical calculations, it could ease the numerical 
effort necessary to perform the tuning, thus reducing the number of required simulations. 

6. Conclusions 



We have presented the results of the non-perturbative tuning of K and zt for the #SF at three 
physical scales and for a range of lattice spacings. This demonstrates that the tuning of these two 
coefficients is indeed feasible, at least in the quenched approximation. Moreover, we observe that 
the tuning of zt and K are nearly independent. Note that even with non-improved Wilson fermions 
in the bulk, K and zt are the only parameters that must be tuned within the #SF setup in order to 
guarantee bulk automatic 0(a) improvement, thus eliminating the need for the bulk counterterm, 
c sw , and for the many operator improvement coefficients necessary in the sSF. 

Our next step is to perform an universality test of this formulation as well as a demonstration 
that automatic 0(a) improvement holds. This can be done by reproducing a variety of quantities 
already computed in the standard setup. A natural candidate would be the computation of the step- 
scaling function of the pseudoscalar renormalisation factor, Zp, which could be compared to the 
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results of [15 1. We recall that the #SF and the sSF are equivalent in the continuum limit, therefore, 
it is not necessary to recompute the entire evolution of an operator. The only quantity that must be 
recomputed is the renormalisation factor at the most non-perturbative scale. 

We also plan to explore whether the value of K^ r determined from the finite volume simulations 
can be used in large volume, preserving the nice scaling behaviour obtained in Refs. Qlq, [l7| ], 
without the need for a large volume determination of K^. A lattice perturbation theory computation 
of d & and Zf is also planned. The final goal is to perform dynamical simulations. 
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